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Q^, Abstract 

-)— > \ 

Ch . We introduce a self-consistent theoretical framework associated with the Schwinger 

^ \ unitary operators whose basic mathematical rules embrace a new uncertainty prin- 

O^' ciple that generalizes and strengthens the Massar-Spindel inequality. Among other 

remarkable virtues, this quantum-algebraic approach exhibits a sound connection 
with the Wiener-Kinchin theorem for signal processing, which permits us to deter- 
J^ ' mine an effective tighter bound that not only imposes a new subtle set of restric- 

l/-^ ■ tions upon the selective process of signals and wavelets bases, but also represents 

CN , an important complement for property testing of unitary operators. Moreover, we 

establish a hierarchy of tighter bounds, which interpolates between the tightest 
^^ . bound and the Massar-Spindel inequality, as well as its respective link with the 

^r^ . discrete Weyl function and tomographic reconstructions of finite quantum states. 

We also show how the Harper Hamiltonian and discrete Fourier operators can be 
combined to construct finite ground states which yield the tightest bound of a given 
finite-dimensional state vector space. Such results touch on some fundamental ques- 
tions inherent to quantum mechanics and their implications in quantum information 
theory. 



1 Introduction 



Initially introduced by Schwinger [1] for treating finite quantum systems cha- 
racterized by discrete degrees of freedom immersed in a finite- dimensional 
complex Hilbert space [2], the unitary operators gained their first immediate 
application in the formal description of Pauli operators. Ever since, an expres- 
sive number of manuscripts [3] proposed similar theoretical frameworks with 
intrinsic mathematical virtues and concrete applications in a wide family of 
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physical systems — here supported by a finite space of states. With regards 
to these state spaces, it is worth mentioning that certain algebraic approaches 
related to quantum representations of finite-dimensional discrete phase spaces 
were constructed from this context in the past [4], and tailored in order to 
properly describe the quasiprobability distribution functions [5] in complete 
analogy with their continuous counterparts [6]. Thus, applications associated 
with the discrete distribution functions covering different topics of particular 
interest in physics — e.g., quantum information theory and quantum com- 
putation [7,8], as well as the qualitative description of spin-tunneling effects 
[9], open quantum systems [10] and magnetic molecules [11], among others — 
emerge from these approaches as a natural extension of an important robust 
mathematical tool. 

Although the efforts in constructing a sound theoretical framework to deal 
with finite-dimensional discrete phase spaces have recently achieved great ad- 
vances (e.g., see Ref. [12]), certain fundamental questions particularly associ- 
ated with the factorization properties of finite spaces [13], uncertainty principle 
[14] and property testing [15] for the unitary operators still remain without 
satisfactory answers in the literature (indeed, some of them represent open 
problems which do not share the same rhythm of progress). In this paper, we 
focus on the problem of deriving a general uncertainty principle for Schwinger 
unitary operators in physics. In what follows, we discuss the relevance of such 
a principle and, subsequently, briefiy review the results obtained by Massar 
and Spindel [14] on this specific subject. 

To begin with, it is necessary to remember that, through an original algebraic 
approach which encompasses the description of finite quantum systems, Weyl 
[16] was the first to describe quantum kinematics as an Abelian group of ray 
rotations in the system space. According to Weyl: "The kinematical structure 
of a physical system is expressed by an irreducible Abelian group of unitary 
ray rotations in system space. The real elements of the algebra of this group 
are the physical quantities of the system; the representation of the abstract 
group by rotations of system space associates with each such quantity a def- 
inite Hermitian form which 'represents' it." With respect to the particular 
case of finite state vector space, one of Weyl's most significant achievements 
was that the observation of pairs of unitary rotation operators obey special 
commutation relations (bringing, as a result, the roots of unity) which are the 
unitary counterparts of the fundamental Heisenberg relations. Moreover, it is 
worth mentioning that such a ray representation of the Abelian group of rota- 
tions can be connected with some representations of the generalized Clifford 
algebra, this fact being thoroughly explored by Ramakrishnan and coworkers 
[17] through extensive studies of certain physical problems. Still within the 
aforementioned Weyl approach for quantum kinematics, let us briefly mention 
that some authors have also addressed the problem of discussing quantum 
mechanics in flnite-dimensional state vector spaces, where the coordinate and 



momentum operators (characterized by discrete spectra) play an essential role 
in this context [18]. 

Although both Weyl and Schwinger's theoretical approaches have made se- 
minal and complementary contributions upon the scope of unitary operators 
in finite physical systems, the relevance of a general uncertainty principle for 
such operators has not been clearly discussed or even mentioned with due em- 
phasis in the past. Reflecting on this, Massar and Spindel [14] have recently 
established a flrst uncertainty principle for the discrete Fourier transform [19] 
whose range of applications in physics covers, among other topics, the Pauli 
operators, the coordinate and momentum operators with flnite discrete spec- 
tra, the modular variables, as well as signal processing. Furthermore, their 
result can also be employed to determine a modifled discrete version of the 
Heisenberg-Kennard-Robertson (HKR) uncertainty principle which resembles 
the generalized uncertainty principle (GUP) in the quantum-gravity frame- 
work [12]. However, if one adopts an essentially pragmatic point of view, cer- 
tain natural questions arise: "Can Massar-Spindel inequality be recognized 
as a 'generalized uncertainty principle' for all flnite quantum states?" If not, 
"What is the reliable starting point for obtaining a realistic description of this 
generalized uncertainty principle?" 

The main goal of this paper is to present a self-consistent theoretical fra- 
mework for the Schwinger unitary operators which embodies, within other 
virtues, an important set of convenient inherent mathematical properties that 
allows us to construct suitable answers for the aforementioned questions. This 
theoretical framework, constituted of numerical and analytical results, can be 
interpreted as a "generalized version" of that one by Massar and Spindel, with 
immediate applications in quantum information theory and quantum compu- 
tation, as well as in foundations of quantum mechanics. Next, we emphasize 
certain essential points of our particular construction process: (i) Numerical 
computations related to a huge number (^ 10^) of randomly generated flnite 
states demonstrate the existence of a nontrivial hierarchical relation among 
the different bounds, the Massar-Spindel inequality being considered in such 
a case as a zeroth-order approximation, (ii) The existence of a tightest bound 
for different dimensions of state vector space leads us to produce a sufficient 
number of formal results related to the Hermitian trigonometric operators (de- 
fined through well-known specific combinations of unitary operators) and their 
corresponding Robertson-Schrodinger (RS) uncertainty principles [20], which 
culminates in the formulation of a new inequality which takes into account the 
quantum correlation effects. This tighter bound represents a new and impor- 
tant paradigm for signal processing with straightforward implications on finite 
quantum states [21] and discrete approaches in GUP [22]. (iii) Numerical and 
analytical approaches [23,24] confirm the special link between the ground state 
inherent to the Harper Hamiltonian and the tightest bound for any Hilbert 
space dimensions. Finally, (iv) the connection with tomographic measurements 



of finite quantum states via discrete Weyl function represents, in this case, a 
tour de force in our investigative journey on unitary operators that allows to 
join both the Weyl and Schwinger quantum-algebraic approaches in an elegant 
way. 

This paper is structured as follows. In Section 2, we fix a preliminary ma- 
thematical background on the Schwinger unitary operators, which allows us 
to discuss the implications and limitations of the Massar-Spindel inequality 
for signal processing. In Section 3, we introduce four Hermitian trigonomet- 
ric operators through effective combinations of the Schwinger unitary ope- 
rators. Together, these operators provide a self-consistent quantum-algebraic 
framework, leading us to determine a new tighter bound for Massar-Spindel 
inequality. Section 4 is dedicated to discuss certain important aspects of the 
tightest bounds and their respective kinematical link with the Harper Ha- 
miltonian. In addition. Section 5 presents an elegant mathematical procedure 
for measuring a particular family of expectation values — here mapped upon 
finite-dimensional discrete phase spaces and related to the unitary operators 
under investigation — via discrete Weyl function. Section 6 contains our sum- 
mary and conclusions. Finally, Appendix A concerns the Harper functions and 
their respective connection with the tightest bounds verified in the numerical 
calculations. 



2 Preliminaries 



In order to make the presentation of this section more clear and self-contained, 
we begin by reviewing some essential mathematical prerequisites related to 
the Schwinger unitary operators. Only then we establish the Massar-Spindel 
inequality and its inherent limitations. 



2. 1 Schwinger unitary operators 

Definition (Schwinger). Let {U, V} be a pair of unitary operators defined 
in a A^-dimensional state vector space, and (Imq), jt'/?)} denote their respective 
orthonormal eigenvectors related by the inner product {ualvp) = -j^oj^^ with 

oj := exp (^). The general properties 

together with the fundamental relations 

U^ = 1, V^ = 1, V«U'' = w^^U^V^ 



constitute an important set of mathematical rules that are related to the 
generalized Clifford algebra [17]. Here, the discrete labels {a, f3, !],(,} obey 
the arithmetic modulo N and (ualvjs) represents a symmetrical finite Fourier 
kernel. A compilation of results and properties associated with U and V can 
be found in Ref. [1]. 



In the following, let p describe a set of physical systems labeled by a finite 
space of states, whereas "^ := 1 — |(U)p and "^ := 1 — |(V)p denote the 
variances related to the respective unitary operators U and V — in this 
case, (U) and (V) represent the mean values of U and V defined in a A^- 
dimensional state vectors space. Since p refers to a normalized density opera- 
tor, the Cauchy-Schwarz inequality allows us to prove that |(U)p and |(V)p 
are restricted to the closed interval [0, 1]; consequently, both the variances are 
trivially bounded by < '^^(v) < 1- Indeed, the upper and lower bounds are 
promptly reached when one considers the localized bases {\ua){ua\}Q<a<N-i 
and {|^^/3)(^^/3|}o</3<Ar-i, i-e-, for a given p = \ua){ua\ ^ % = ^ and % = \\ 
otherwise, if p = |u^)(t;^| ^ Xj = 1 and X' = 0. Furthermore, note that "V^ 
and "V^j are invariant under phase transformations, namely, U — )■ e^'^U and 
V -^ e'^y for any {cp, 6} G M. 



2.2 Massar-Spindel inequality 



This inequality is based on the Wiener-Kinchin theorem for signal processing 
and provides a constraint between the values of (V ) (correlation function) 
and (U ) (discrete Fourier transform of the intensity time series). According 
to Massar and Spindel: 'This kind of constraint should prove useful in signal 
processing, as it constrains what kinds of signals are possible, or what kind of 
wavelet bases one can construct.' We state this result in the theorem below 
(proved in the supplementary material from Ref. [14]), for then proceeding 
with a numerical study of its content and first implications. 

Theorem (Massar and Spindel). Let U and V denote two unitary operators 
such that U V = e VU and U V = e^ VU with <1> G [0,'7r). The variances 
y^j and y— — here defined for a given quantum state p and limited to the closed 
interval [0, 1] — satisfy the inequality 



\ 



(1 + 2A)i%r- >A\i-Yjj- i\ 



v7 



where A = tan (-j) . The saturation is reached for localized bases. 




This theorem leads us, in principle, to consider the different possibilities of 
connections between the Schwinger unitary operators {U, V} (defined in the 



previous subsection) and {U,V}. A first immediate link yields the relation 
U = U and V = V for $ = — ^, which implies in the apparent violation 
of Eq. (1) since $ ^ [0,7r). The second connection establishes the alternative 
relation U = V and V = U with $ = ^, this result being responsible for 
validating the Massar-Spindel inequality. It is important to mention that the 
apparent problem detected in the first situation can be properly circumvented 
making $ — ?■ — $ and A — )■ —A with $ = ^ fixed. That is, the modified 
bound 

(1 - 2A)%ji<, > A\l -%- %) 

holds for any A = tan ( -^ j and integer A^ G [2, oo). 

For the sake of simplicity and convenience, let us now introduce the shift 
operator AO = O — (O) with (O) ^ for a given arbitrary unitary operator 
O and density operator p, which leads to define ^ its non-unitary counterpart 
as follows [25]: 

50 := —— = — ^— ^. (2) 

Following, it is straighforward to show that both the variances i^a and i^so S'-re 
related through the expressions 

rso = -^ (0<r,,<oo) and r, = ^^, (3) 

which substituted into inequality (1) for U = V and V = U gives 

ysvy5v>e^ with e = ---. (4) 

1 -\- A 



To illustrate Eq. (4) and corroborate the analytic results obtained by Massar 
and Spindel, let us introduce the parameters Ssv = ^'^^v and Sgy = e^^Ysv 
such that Ssu<Ssv > 1- This particular inequality defines a region in the two- 
dimensional space limited by the rectangular hyperbola Ssv = S^^ that pre- 
serves the original equation l^u^v > ^^- Figure 1 shows the plots of Ssy 
versus Ssv for approximately 10^ states of randomly generated p, within the 
visualization window < Ssv,Ssy < 10, with (a) A^ = 2 (e — )■ 1), (b) A^ = 3, 
(c) A^ = 4, and (d) A^ = 5. Note that, excepting picture (a), all the subsequent 



^ It is important to note a certain level of arbitrariness in this definition because 
there is no especification whatsoever of which density operator p is used to evaluate 
the mean value in the denominator of Eq. (2). Henceforth, this arbitrariness will 
be removed by exploring the fact that, in our computations, the operator 60 will 
always appear as an argument of a variance function — in such a situation, we will 
understand that the expectation value (O) is computed with respect to the same 
state used in the computation of {60). 
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Fig. 1. Plots of Ssv versus ^^u for different values of A^ — 2 (a), 3 (b), 4 (c), 
and 5 (d) — with approximately 10^ randomly generated states, illustrating the 
Massar-Spindel inequality (^^u^^v ^ 1) for the Schwinger unitary operators. Note 
that < Ssu,<Ssv < 10 characterizes a particular visualization window which leads 
us to guess the existence of a tightest bound (see dashed curves) for each value of 
dimension N, since the distance between the saturation curve SsuSsv = 1 (solid 
line) and the cloud of points — generated by numerical calculations — increases for 
N > 2. The dot-dashed curves showed in pictures (b,c,d) represent the intermediate 
inequality ^^u^^v > 1+sin (^) , this result being considered as a first approximation 
to our initial intents. 



cases exhibit a gap between the distribution of states and the rectangular hy- 
perbola, the size of such a gap being dependent on the value of A^ (such an 
evidence motivates the search for a tightest bound that corroborates the nu- 
merical calculations). In fact, the dashed lines in (b,c,d) describe hyperbolic 
curves of the form SsvSsv = x, where the value of x was chosen as the smallest 
value of Ssv^sv amongst the ones computed with the randomly sampled states. 
Later in this paper, a more rigorous procedure for obtaining such values will 
be outlined. Finally, let us briefly mention that the dot-dashed lines exhibited 
in these pictures correspond to the intermediate result SsvSs^ > 1 + sin f ^j, 
which will be properly demonstrated in the next section. 



3 A Hierarchy of Tighter Bounds 



In the first part of this paper, we established a basic theoretical framework re- 
lated to the Schwinger unitary operators where the Massar-Spindel inequality 
and its inherent limitations occupied an important place in our discussion on 
uncertainty principles for physical systems labeled by a finite space of states. 
At this moment, let us clarify some fundamental points raised by those results: 
(i) the aforementioned state spaces consist of A-dimensional Hilbert spaces; 
(ii) the Massar-Spindel inequality represents a "zeroth-order approximation" 
in the hierarchy of uncertainty principles; and finally, (iii) the results obtained 
from the numerical calculations reveal certain unexplored intrinsic properties 
of some finite quantum states [26,27] with potential applications in quantum 
information theory and quantum computation. In this second part, we begin 
the construction of a solid algebraic framework based on the RS uncertainty 
principle, which leads us, in a first moment, to determine a self-consistent set 
of results for the unitary operators U and V that permits to generalize the 
Massar-Spindel inequality. Indeed, these results represent an important tool 
in our search for tighter bounds (see numerical evidence exhibited in Fig. 1), 
whose intermediate uncertainty principles will constitute a hierarchical rela- 
tion between the Massar-Spindel inequality and the tightest bound. 



3.1 Quantum- algebraic framework 



Let {A, B} denote a pair of Hermitian operators defined in a A^-dimensional 
state vectors space which obey the RS uncertainty principle [20] 

r.rB>^^(A,B) + l|([A,B])|% (5) 

where "^(A, B) := (|{A, B}) — (A)(B) represents the covariance function, and 
[A, B] ({A, B}) corresponds to the commutator (anticommutator) between A 
and B. ^ Next, since U and V are generally non-Hermitian operators, let us 
consider the cartesian decomposition of an arbitrary non-Hermitian operator 



^ It is important to emphasize that, according to Cauchy-Schwarz inequahty, the 
covariance function ^(A, B) is restricted to the closed interval [— V^^a^, \/^a^] , 
namely, |^(A,B)| < v^a^- Indeed, for a given operator C := A — "-^ — ^B with 

1^/0 and '^ > 0, it turns immediate to obtain the relation i^c = i^A y ' , 

which demonstrates the previous result for "^(AjB). Moreover, let us also state a 
very useful result for both the commutation and anticommutation relations between 
A and B, that is |([A,B])|2 + |({A,B})|2 = 4|(AB)|2. 



O into its 'real' and 'imaginary' parts as follows [28]: 

Note that Re[C?] and Im[C?] represent two Hermitian operators that comply 
with the RS uncertainty principle and allow to introduce, in particular, the 
cosine and sine operators through the Schwinger unitary operators. 



Definition. Let {Cu, Su, Cy, Sy} denote four Hermitian operators written in 
terms of simple combinations associated with U and V, i.e., Cu := Re[U], 
Su := Im[U], Cy := Re[V], and Sy := Im[V]. The commutation relations 
involving these operators exhibit a direct connection with certain anticommu- 
tation relations: 

[Cu, CyJ = lA |bu, Syj , [Cu, SyJ = — 1/1 |bu, Cyj , 

[Su, SyJ = m jCu, Cyj , [Su, CyJ = —lA jCu, Syj , 

where the parameter A was previously defined in Eq. (1). These results lead us, 
in principle, to conclude that partial information on two particular elements of 
the set is not complete, since the complementary elements are also necessary 
to fully characterize the commutator [U, VJ. In this sense, let us now consider 
the four RS uncertainty principles below: 



rCurCy><^^(Cu,Cy) + i|([Cu,CyJ)|% (6) 

^Cu rSy > <^2(Cu, Sy) + ^ I ([Cu, SyJ) I' , (7) 

rSu^Cy > <^^(Su, Cy) + i I ([Su, Cyj) I' , (8) 

^Su^Sy > <^2(Su, Sy) + i I ([Su, SyJ) 1^ . (9) 

In addition, the extra result CL^n + SL^n = 1 resembles a well-known mathe- 
matical property associated with the trigonometric functions cosine and sine. 
For this reason, these Hermitian operators will be henceforth termed 'cosine' 
and 'sine' operators, whose respective variances can be shown to obey the 
mathematical identity 

(^Cu+^Su)(^Cy + ^Sy) = >;>;• 

To make complete this definition, it is interesting to observe that certain com- 
binations of (Cu(y))^ and (Su(y))^ also yield the additional result 

((Cu)^ + (Su)^) ((Cy)2 + (Sy)2) = (1 - %){! - %), 



which proves itself useful in our subsequent calculations. 



Next, by means of mathematical remarks, we establish an important set of 
other useful results for the sine and cosine operators, whose relevance is in- 
trinsically connected with the hierarchy relations involving the uncertainty 
principles related to U and V which generalize the Massar-Spindel inequality. 
It is worth mentioning that some proofs demand a logical sequence of algebraic 
manipulations to be detailed in the text. 



Remark 1. Although the sine and cosine operators are genuinely defined 
for unitary operators, we shall employ, in due time, the same definition for a 
general operator X as well. In this general case, it can be easily demonstrated 
that the following properties hold: (C|) + (S|) = (|{X,Xt}), (Cx)=^ + (Sx)^ = 
I (X) p, and fcx + '^Sx = ^(X, X^"). Note that for a normal operator N (which 
satisfies [N, N^] = 0), the covariance function between N and N^ matches the 
variance of N, namely, Yq^ + Ys^ = %i- 



Remark 2. Let us initially consider the sum of all aforementioned RS uncer- 
tainty principles for the sine and cosine operators, as well as the connection 
between the commutation and anticommuation relations of such Hermitian 
operators. Adequate algebraic manipulations allow, in principle, to obtain an 
inequality for the product YuY/ with remarkable mathematical features, i.e.. 



1 + A^)%% > (1 + A2)J^(U, V) + A^ nilJ, V) (10) 



where 



J-(U, V) = ^2(c„, Cv) + ^'(Cu, Sv) + ^'(Su, Cv) + ^'(Su, Sv), 

n{v,v) = |(CuCv)|' + |(CuSv)|' + |(SuCv)|' + |(SuSv)|' . 

Note that J-" and Ti show a nontrivial dependence on the unitary operators U 
and V, which will be properly discussed in the subsequent remarks. Moreover, 
if compared with Massar-Spindel inequality (1), such a result yields subtle 
additional corrections that will depend explicitly on the dimension of the state 
space. 



For completeness sake, let us now establish a first numerical evaluation on the 
functions J-" and "H. Figure 2 exhibits the plots of J-' versus "H for (a) A^ = 4 
and (b) A^ = 5 with the same number of normalized random states used in the 
previous figure. Since the solid line represents J-" = "H in both the situations, 
this preliminary numerical search demonstrates a greater contribution coming 
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(b) 




0.2 ^ 0.3 




0.2 ^ 0.3 



Fig. 2. Plots of J^(U, V) versus ^(U, V) with 10^ states of randomly generated p 
for two different values of dimension iV, namely, (a) A^ = 4 and (b) N = 5. The solid 
line corresponds to the case J- = Ti in both the pictures, where the limit situation 
J- = Ti = describes the localized bases. Our numerical computations demonstrate 
that, in principle, the contributions due to Ti are more significant than those coming 
from J^ for a large number of random states. 

from "H than from J-" for most states of low dimension A^. This fact appeals for 
a detailed formal investigation on the origins of such contributions, in which 
each term of J-' and Ti would be examined separately. 



Remark 3. Important mathematical properties of 'H(U,V) are easily at- 
tained examining each of its terms separately: for instance, to calculate the 
modulus squared of (CuCv), we initially expand the cosine operators Cu and 
Cv in terms of the Schwinger unitary operators U and V; the next step consists 
in decomposing the resulting terms into their real and imaginary Hermitian 
parts, whose final expression assumes the form 



(C„Cv) = - ((C„v) + (C„vt) + (C„tv) + (C„tvt)) 

+ ^((S„v) + (S„vt) + (S„tv) + (S„tvt)). 

Once the remaining terms are obtained in the same way, some algebraic sim- 
plification gives rise to the following expression for "H: 



n{v, V) = 1 ((c„v)' + (c,vt)' + (c„tv)' + (c.tvt)') 

+ I ((Sw)' + (S„vt)' + (S„tv)' + (S„tvt)' 

However, this result does not represent a convenient form of "H since the proper 
summation of (Cq)^ + (So)^ = l — Yo with (C^) + (S^) = 1 (in this situation, 
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one considers all possible matches of O = UV, UV''', U^V, U^V''') permits us 
to derive a simplified expression for such a function that includes the variances 
^v and l^vt- Indeed, since "^tv = "^vt and "^tvt = '^v, it turns immediate 
to prove that 

?/(u, V) = 1 - 1 (r„, + r„vt) = |(uv)p + |(uvt)|2 ^^^^ 

which is invariant under the transformations U — )■ e"^U and V — )■ e'^V for 
{ip,e} eM. 



Remark 4. Note that Eq. (11) represents a lower bound for any element of ^ 

{|(c„Cv)lM(c„Sv)lM(s„Cv)lM(s„Sv)p} > ^^(u, v), 

this fact being discussed by Massar and Spindel [14] through different mathe- 
matical arguments. In fact, the authors demonstrated that for a given choice 
of phases of the Schwinger unitary operators, the restrictions (Cu(v)) ^ ^+ =^ 
(Su(v)) = select a set {p} of finite quantum states for which the inequality 



|(CuCv)| > ^{i-%){i-%)-^/f^ 

can be formally verified and also numerically tested. Despite the correlations 
between U and V do not appear in such an expression, the comparison with 
"H is unavoidable in this case, since correlations represent an important quan- 
tum effect that deserve our attention. The saturation is reached in both the 
situations for localized bases. 



Remark 5. Let us now decompose J-' into three terms J-'i, T2 and J-3, whose 
different contributions will be formally calculated with the help of the mathe- 
matical procedure sketched in the previous remarks. We initially consider the 
term 

7^ = \ (({C„ Cv})^ + ({C„ Sv})^ + ({S„ Cv})^ + ({S„, Sv})^ 



^ Let ||U||j^g = ||V||jjg = \fN and ||Cu(v)|Lo = ||Su(v)||ttc = a/y characterize the 
Hilbert-Schmidt norms associated with the Schwinger unitary operators and their 
respective related trigonometric operators, where ||X||jjg := y^Tr [XtX] defines the 
aforementioned norm [28]. The further mathematical property 



illP P II IIP c; II ||c p II ||c c II \ < 

\ II'-^u'-^vIIhs ' II^u'-'vIIhs ' II'-'u'--^v|Ihs ' II'^u'-'vIIhs | — 



N 

T 



represents an effective gain in this stage, since it brings out relevant information on 
the different products of cosine and sine operators used in the text. 
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responsible for contributions associated with squared mean values of all anti- 
commutation relations involved in T through the covariance functions. So, let 
({Cu, Cy}) admit the form 



({C„ Cv}) = ^ [(1 + ^) ((C„v) + (C„tvt)) + (1 + a;*) ((C„vt) + (C„tv))] 
+^ [(1 + u) ((S„v) + (S„tvt)) + (1 + c^*) ((S„vt) + (S„tv))] 

in complete analogy with (CuCv)- Writing the remaining mean values in an 
analogous way, subsequent computations allow to prove that T\ is connected 
with ?^(U, V) through the relation J=x = cos^ (f ) ?{ for $ = ^ fixed. 

Following, in what concerns the second term 



J-2 = (Cu) (Cv) ({Cu, Cv}) + (Cu) (Sv) ({Cu, Sv}) 

+ (Su)(Cv)({Su, Cv}) + (Su)(Sv)({Su, Sv}), 

it can be expressed by means of the convenient form 



2J-2 = Re[(U)(V)]Re[(l + a;)(UV)] + Re[(U)(Vt)]Re[(l + uj*)(UY^)] 
+Im[(U)(V)]Im[(l + a;)(UV)] + Im[(U)(Vf)]Im[(l + uj*){UY^)] 

which represents an important formal result in our calculations. 

Finally, let us briefly mention that J-3 coincides with (1 — '^)(1 — "K) since 
it is equivalent to the product ((Cu)^ + (Su)^) ((Cv)^ + (Sv)^); consequently, 
the function 

^(U, V) = J-i(U, V) - J-2(U, V) + J-3(U, V) (12) 

can be immediately obtained. Note that the nontrivial dependence of Eq. (10) 
on the Schwinger unitary operators is completely justified in these remarks. 



This set of mathematical remarks establishes a first solid algebraic framework 
for the unitary operators U and V, whose intrinsic virtues lead us to formulate 
a theorem which generalizes the Massar-Spindel inequality (1) in order to 
include the quantum correlation effects between the aforementioned operators. 
In fact, this theorem consists of an initial compilation of efforts in our future 
search for the tightest bound, where the correlation function has occupied an 
important place in the investigative process. 
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■ I 

Theorem. Let U and V be two unitary operators defined in a N -dimensional 

state vector space which satisfy the commutation relations [U, V] = (1 — u;)UV 
1 and [U, Vt] = (1 -w*)UVt for to := exp (^). Moreover, let ^u := 1 - |(U)|2 

and Yy := 1 — |(V)p represent the respective variances such that < '^(v) ^ 1- 

The inequality 
\ %% > -F(U, V) + sin^ (^) mV, V) (13) 

yields a new bound for 'f^'Vy, where the quantum correlation effects related to 
the unitary operators are taken into account. Note that T and % were pre- 
cisely defined and formally studied in the previous mathematical remarks, the 
saturation T = Ti = being attained for localized bases. 

Next, let us determine some additional results focused on the Hermitian ope- 
rators {C^u! S^U) C^v, S^v} in order to yield a set of specific intermediate ine- 
qualities whose hierarchical relations correspond to a solid bridge towards the 
tightest bound. 



3.2 Hierarchical relations 



We start this subsection stating a first important result for the sine and cosine 
operators {C^u; S^u, C^v; S^v} defined, in turn, in terms of the non-unitary 
counterparts 6\J and 6\. This particular result shows how the mean values of 
their commutation and anticommutation relations are linked with determined 
correlation functions, namely, 

{[Cs^,Csy]) = iA{{Ss^,Ss.}) = 2iA^{Ss^,Ss.), 

([S,u,C5v]) = -iA({C,u,S5v}) = -2iA^{Csu,Ss.), 

([S,u, S^v]) = iA (({C^u, C^v}) + 2) = 2iA [^{Cs^, Cs.) + 1] . 

Following, let us complete this set of results with relations that connect all 
the covariance functions previously defined with the original framework: 

J-3 ^(C^u, C,v) ='^(Cu, Cv)(Cu)(Cv) + ^(Su, Sv)(Su)(Sv) 

+^(Cu, Sv)(Cu)(Sv) + ^(Su, Cv)(Su)(Cv), 

T3 ^(C^u, S^v) = — ^(Cu, Cv) (Cu) (Sv) + "^(Su, Sv) (Su) (Cv) 

+^(C„, Sv) (C,) (Cv) - <^(S„, Cv) (S„) (Sv) , 

J-3 ^(S^u, C5v) = -^(Cu, Cv)(Su)(Cv) + <^(Su, Sv)(Cu)(Sv) 

-^(Cu, Sv)(Su)(Sv) + <^(Su, Cv)(Cu)(Cv), 
J^3^{Ss^, S,v) =^(Cu, Cv)(Su)(Sv) + ^(Su, Sv)(Cu)(Cv) 

-^(Cu,Sv)(Su)(Cv)-^(Su,Cv)(Cu)(Sv). 



14 



These equations enable us to rewrite T in the compact form 

-— =^ (C^u, C^v) + ^ (C^u, S^v) + ^ (S^U) C^v) + ^ (S^U! S^v) 
-^3 



(14) 



while (l^^u + '>^iu)(^<sv + ^4v) ~ ^u'^^v brings out a completely analogous 
expression to that established for YuYy] consequently, 



r.A^>^ + si„=(^) 



7r\ 2i_ 
•F3 



(15) 



represents an alternative form of Eq. (13) since Y^jY-m = T'iVsv"^^ 



Sv 



Now, let us rewrite Eq. (15) according to definitions employed for Ssu and Ssv 
in the previous section. 



'-'<5u'-'<5v ^ 



, 27r 
l + sin(- 



CSC 



n 



7 n 



N) J-3 J-3 



(16) 



From this result, it is relatively easy to derive a number of simpler (but looser) 
inequalities: for instance, given the non- negativity of ^ [cf. Eq. (14)], the first 
term in the second bracket can be dropped to give 



'-'5u'-'<5v ^ 



r /27r\l 


n 


1 + sin — - 




[ \N )\ 


•Fa' 



which corresponds, in such a case, to the HKR uncertainty principle; alterna- 
tively, -^ can i 
order to yield 



tively, -^ can also be dropped (for the same reason, i.e., its non-negativity) in 



Consequently, such inequahti 


1 + sm -- 

[ \nJ\ 

es are then co 


CSC -- -—. 
mbined in order to give a tighter 


one. 


1 + s 


/27r\1 


max 


" 2 f^\ J^ ^' 


(17) 



At this point, one is led to ask if there exists a general ordering relation between 
csc^ ( -^ j ^ and ^. To answer this question, let us establish a set of numerical 
calculations associated with the different values of dimension A^, where, for 
each specific case, one has approximately 10^ randomly generated states. In 
order to make the presentation of these numerical results more self-contained. 
Fig. 3 depicts the plots of csc^ f -^J ^ versus ^ for A^ = 2, 3, 4, 6 and 20. In this 
case, note that the maximum can indeed arise from either terms, depending 
on the particular chosen state p. As a rule of thumb, one has the following: 
for low dimensional states (e.g., A^ = 2,3), the term y- usually dominates; 

however, as the dimension A^ increases (A^ > 6), csc^ (j^]y- becomes the 
usually dominant term. Finally, if one considers A^ = 4, it is visible that there 
is not a usually dominant term — in this case, the result of the optimization 
is strongly dependent on the particular input state. 
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100 




100 



Fig. 3. Plots of csc^ (■^) y- versus ^ with, at least, 10^ randomly generated states 
for each different value of dimension N. Thus, yellow points depict the N = 2 case, 
while red, grey, blue, and green points describe, respectively, the A^ = 3, 4, 6, and 20 
situations. Note that (0, 0) is associated with localized bases in this picture, whereas 
the solid line corresponds to the particular case csc^ (;^) ^ = ^. It is interesting 
to observe how csc^ (-^) -i- and -^ behave when different dimensions of state vector 
space are considered in our numerical computations. 

As it turns out, inequality (17) can be further strengthened by the addition of 
yet another (less trivial) lower bound of csc^ f -^ ) y- + y- to the maximization 
at hand. Just as -^ could be written as in Eq. (14), we can also show that ^ 



admits the form 



^3 



n 2 f'^ 

— = sec — 



l + 2<^(C5u,C5v) + ^ 



With that in mind, some straightforward manipulation yields 



7r\ J^ H 



sec 



TC 



N 



— l+Xr + CSC^ TT he 



--='!) (^ 



X 



TC 



N 



>sec2(^)(l + a;)2 + csc2(^^)a;2 



for X = ^{Csv, Csv), where the inequality follows from fact that y- > x"^ [cf. 
Eq. (14)]. A mathematical virtue of this lower bound is that it depends only 
on the covariance ^{Csv, C<5v), whereas ^ and ^ depend on the covariances 
of all combinations among C<5u, Csv, Ssu and S^v As a result, we may now 
write 



'-'(5u'-'<5v ^ 



1 + sin 



27r 

iV 



max 



I + g{x;N), csc^ 



vr 



J" U 



Nj :f.' :f. 



(19) 
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with g{x; N) := 4csc2 (^) x2 + 2sec2 (^) x + tan^ (^). Note that l + g{x; N) 

coincides, in such a situation, with sec^ f-^j (l + a;)^ + csc^ f -^j x^. Once again, 
numerical calculations demonstrate that there is not a general ordering be- 
tween the arguments of the maximization of such an equation. 

It is worth stressing that both bounds of Eqs. (16) and (19) explicitly depend 
on the Hilbert space dimension, as well as on the particular state under con- 
sideration. Henceforth, we further loose the bound (19) to provide yet another 
bound, but now a state-independent one (i.e., solely dependent on the Hilbert 
space dimension). In this way, let us initially consider the inequality 



^Sv^Sy ^ 



, 27r 
l + sin(- 



[l + gix;N)] 



So, for a given A^, Q{x; N) describes a parabola with upwards concavity and 
minimum value equal to (y-coordinate of the vertex), namely, g{x; N) > 0, 
which implies that 

SsvSsv > 1 + sin C^) . (20) 



NJ 

Despite of all the mathematical assumptions used in the relaxation process 
of inequality (16), the bound above is still tighter than the one proposed in 
Ref. [14], which is now trivially proved by disregarding the sine function in 
the inequality (20), that is, Ss^jSsy > 1. 



alem | 



Hierarchy. The tightest bound of the product SsySsy is particularly relate' 
an underlying minimization problem — see Appendix A. Although this proble: 
seems simple at first glance [23] , let us now establish a set of inequalities which 
characterizes an important hierarchical relation among certain formal results 
mentioned in the body of the text: 



I SsvSsv > IHi > IH2 > 9^3 > 9^4 > 1, (21) 

where 9^i = tightest bound, 9^2 = RHS of Eq. (16), 9^3 = RHS of Eq. (19), 
and finally, 9^4 = RHS of Eq. (20). Note that ^^u^^v ^ 1 describes the loosest 
bound, which is precisely the inequality of Massar and Spindel [14]. 



4 The tightest bound 



Massar and Spindel [14] relied on the works of Jackiw and Opatrny [23] to 
provide a numerical recipe that implicitly produces the tightest bound for the 
cloud of points depicted in Fig. 1. In this section, we review such a procedure 
while applying it to explicitly construct closed-form expressions for the tightest 
bound here related to physical systems characterized by a low-dimensional 
space of states. 
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Generally, the tightest bound can be obtained by the following reasoning: 






.-rr 



(0)^ 



'^^u 



In such a case, the super-index (0) indicates that the corresponding normalized 
discrete wavefunction is associated with the nondegenerate ground state \ipo), 
which minimizes the product "^1^. In fact, this ground state will also minimize 
%u'^5v, since Yso is monotonically increasing with Yo [cf- Eq. (3)] — this 
observation justifies the second equality. Besides, both the variances evaluated 
with respect to the ground state satisfy "^^ = l^v , which justifies the third 
equality. 



Now, let us establish a sequence of mathematical steps based on the respective 
eigenvalues and eigenvectors of the Harper Hamiltonian [24] 



H = -sin(^)Cu-cos(^)C, 



for < ^ < |. Initially, we look for the smallest eigenvalue of such a Her- 
mitian operator, as well as for its corresponding eigenvector in a given fixed 
A^-dimensional state vector space. Since the eigenvalues and eigenvectors re- 
lated to the Harper Hamiltonian are dependent on the angle 6, let us consider 
that eigenvector evaluated at the first step in order to estimate the maximum 



- 



- which exactly coincides with the 
ar, such a mathematical procedure 



of cos(^) I (U) I +sin(^) I (V) | for all 9 e 

smallest eigenvalue in this case. In particu 

allows us to obtain, through numerical evaluations, the value oi 9 = ^ for 

any dimension A^, and also to characterize the ground state \iIjq) = \0)n with 

well-defined mathematical properties (see Appendix A for possible connection 

with Harper functions). 

Table 1 illustrates the hierarchical relation depicted in Eq. (21) where, in par- 
ticular, certain numerical results directly related to the analytical calculations 
performed for {9^1,9^2,5^3,5^4} with A^ G [2,6] are exhibited — closed-form 
expressions for S^J can be viewed in Appendix A and their respective squared 

(0)2 

values 9^1 = S^^ compared with those results previously obtained in Fig. 1. 
Note that for higher dimensions, some numerical methods (for instance, La- 
guerre method or even Newton-Raphson method) should be applied in order 
to obtain approximate numerical values for S^J (this statement is supported 
by the Abel-Galois irreducibility theorem [29], which asserts that polynomial 
equations of degree > 5 do not produce, in general, algebraic solutions, only 
numerical solutions). 



Table 1 

Numerical values for {9^j}i<i<4 considering the ground state {|0)Ar}2<Ar<6 described 
in Appendix A. It is important to stress that such results were inferred from their 
respective exact algebraic counterparts, which illustrate, in principle, not only the 
hierarchical relation depicted in Eq. (21) but also the quantum-algebraic framework 
developed in the previous sections. In addition, note that eigenvalues extracted from 
the Harper Hamiltonian for A^ > 6 do not yield easy-to-compute expressions for the 
tightest bounds 9^i, this fact being supported by Abel-Galois irreducibility theorem 
for polynomial equations. 

N mi m2 iHs 5^4 



2 1 1 1 1 

3 ^ 3.254 Rs 2.182 ^ 1.895 ^ 1.866 

4 4 3 2 2 

5 « 3.781 « 3.469 ?« 2.987 « 1.951 

6 ^ 3.348 Ri 1.915 ?« 1.915 ?« 1.866 



5 The connection with discrete Weyl function 



How the discrete Weyl function [5] can be employed to measure a particular 
family of expectation values — for instance, (U"V^) for {a, (3} G Z^ — here 
mapped into finite-dimensional discrete phase spaces? To answer this specific 
question, we initially recall certain basic mathematical tools which correspond 
to the central core of that theoretical formulation presented in Ref. [12]. This 
procedure will lead us to establish a parallel quantum-algebraic framework for 
those results obtained in Section 3, whose connection with the discrete Weyl 
function represents a first step towards effective experimental measurements 
via tomographic reconstructions of finite quantum states [8] . Throughout this 
section, we will assume A^ odd. ^ 

The particular set of A^^ operators { A(/i, i')}fj.,u=~e,...,e characterizes a complete 
orthonormal unitary operator basis which leads us to construct all possible 
kinematical and/or dynamical quantities belonging to a given A^-dimensional 
state vector space. For instance, the decomposition of any linear operator O 
in this basis assumes the expression 

= ^ E ^(^,^)A(/i,z/), (22) 

where (^(/i, z/) = Tr[A(yU, z/)0] represent coefficients evaluated through trace 



^ For completeness reasons, it is important to stress that even dimensionalities can 
also be dealt with simply by working on non-symmetrized intervals. 
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operation and 



A{^^,ly)■.= ^ E ^'^""-^'^^{V^O (23) 



defines the aforementioned operator basis liere expressed in terms of a discrete 
Fourier transform of the displacement operator 



7=-£ 7=-^ 



In such a case, note that to^^"^ ^^^ consists of a specific phase whose argument 
satisfies the mathematical rule 2{2~^riC,} = f]C '^ kN for all k G Zjv; further- 
more, note that these labels assume integer values in the symmetric interval 
[-i, i] with i=^ fixed. 

According to expansion (22), the decomposition of any density operator p in 
the mod(A^)-invariant unitary operator basis (23) has as coefficients the dis- 
crete Wigner function #^(/x, z/) := Tr[A(/i, z/)p], which leads us, in principle, 
to establish an analytical expression for the mean value (O) = Tr[Op], that 
is 

(0) = 1 E ^(/i,^)^p(^,^). (24) 

In what concerns to Wp{fi, u), it is particularly worth mentioning that such a 
function is connected to the Weyl function y^pij],^) := Tr[D(r7,,^)p] by means 
of a mere discrete Fourier transform, and its complexity basically depends on 
the initial quantum state adopted for the physical system under investigation. 

To illustrate the mathematical steps used in the evaluation of ^(/i, z/), let us 
consider those specific combinations of cosine and sine operators exhibited in 
Section 3, as well as the intermediate result 

Dt(r7,OU"V'^D(r7,0 = w^^+^^U^V^ 

since the trace operation 

will be necessary in the next steps. Thus, after some repetitive calculations, 
we achieve the set of formal expressions 
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N 



Tr[D(,, OC„Cv] = ^ c-i-^'^^i (C + Clo) (--4?i,o + -"€i 



N 



Tr[D(,, OC„S,] = i ^ c-i-^'^^i (4^1,0 + €io) (--4?i,o " -^4"l 



N 



Tr[D(r,, OS„C.] = -i ^ c.-i-^^'^> {5[% - 4^' ,o) (c^-4?i,o + -''4''i,o 



4 

which represents, in such a case, the dual counterpart of ^(/i, u) — the super- 
script [A^] on the Kronecker delta denotes that this function is different from 
zero when its discrete labels are congruent modulo A^. 

The second and last step consists in computing, for each case above, its respec- 
tive discrete Fourier transform, which basically depends on the intermediate 
result 



Tr[A(/i,z/)U"V 



a\7-l3] _ , B^J.+au-a|3+{2-''a|3} 



U) 



So, the final expressions admit the simple forms 



(SuCv) (/u, z/) 



1 

4 
i 

~4 
i r 



(CuCv) (/i, i^) = - lico^-" + w^^+'j tu'-^^^^ + {cu>'+'' + CO-''-') w^i+^^-n 
(CuSv) (/i, ly 



^M-- _ ^-f^+A c^l-{2-^} 



^M+-_^-M-^\ -l+{2-i} 



^^_, _ ^_^+A ^l-{2-} _ /,+u _ ^-^-u\ ^-l+{2-} 



(SuSv) (/i, iy) = j [{oo'^-'' + w"^+'^) 0;^-^=^"'^ - (a;^+" + o;"^-") w-^+^^-n 



where (AB) (/i, z/) = Tr [A(/i, z/)AB] denotes the above mapped expressions. 
Consequently, with the aid of Eq. (24), the associated mean values can be 
promptly obtained: for instance, if one considers (U"V^), we easily reach the 
important result 

(U"V^) = c^-^'^+^^-ia/si^^^^ _py (25) 

Now, let us pay attention to (CuCv) and its particular hnk with the discrete 
Weyl function. 



\CuCv) 



y-l,-l) + ;rp(l,ljja; 
^(1, -1) + ^(-1,1)1 cu- 



l-{2-i} 



-l+{2-l} 



(26) 



This remarkable result states that, for a given p and its respective discrete 
Weyl function '^p{ri,^), (CuCy) can be promptly inferred from tomography 
measures upon such a function at specific points of the dual A^-dimensional dis- 
crete phase space. Besides, if one applies the triangle inequality for |(CuCv)|, 
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both the upper and lower bounds can be easily established in this context, ^ 

|(c,Cv)| I \ (|^(-i, -1) + ^(1, 1)1 ± |#;(i, -1) + ^(-1, 1) 

Similar procedure leads us to prove that | (SuSy) | shares exactly the same 
bounds, while the remaining quantities are bounded by the inequalities 



<1 

4 



{|(C„Sv)|, |(S„Cv)|} I - [\W,{-1, -1) - W,{1, 1) 



\{\W,{l,-l)-W,{-l,l)\). 



± 

4 

As a last note, let us mention that 7^ and Yy can also be expressed in terms 
of specific WeyHunctions with the help of Eq. (25), i.e., % = l- |#^(l,0)p 
and rv = 1-1^(0,-1) p. 

In particular, it is worth stressing that the compilation of results obtained in 
this section allows us to rewrite inequality (13) into a new quantum-algebraic 
framework, since J-'(U,V) and 'H(U,V) now depend on specific combina- 
tions of discrete Weyl functions. From an experimental point of view, this 
observation represents an effective gain towards tomographic measurements 
in iV-dimensional discrete phase spaces of such a generalized inequality for 
any finite quantum state p [8]. 



6 Concluding remarks 



Through a well succeeded concatenation of efforts in a recent past [5,12], we 
have made great advances (from a theoretical point of view) in constructing 
certain sound quantum-algebraic frameworks for finite-dimensional discrete 
phase spaces. Since unitary operators represent the basic constituent blocks 
of these theoretical frameworks (in particular, the Schwinger's approach for 



^ Let A and B characterize two general matrices of same size, as well as p denote 
the density matrix related to a physical system described by a finite-dimensional 
state vector space. The additional inequality |Tr[ABp]p < Tr[AATp]Tr[BBTp] (see 
Ref. [27, page 230]) establishes a new upper bound for Eq. (26) since |(CuCv)|^ < 
(C2)(C2), where 



1 1 



(C^) = - + - #;,(2,0) + ^p(-2,0) 



and 



2 4 



(C^) = ^ + 7 #;,(0,2) + ^,(0, 



2 4 

This complementary result basically represents a step forward in our comprehension 
on hierarchical relations involving those means values listed in Remark 4. 
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unitary operators [1]), let us focus our attention on some real and effective 
gains obtained from this paper whose formal implications deserve to be care- 
fully discussed. 

• The Massar-Spindel inequality does not explain the tightest bounds exhibi- 
ted in Fig. 1 for the product SsijSsv when dimensions A^ > 3 are considered 
in the numerical evaluations. This concrete evidence suggests the implemen- 
tation of an effective search for different inequalities and new finite quantum 
states whose implications lead us not only to obtain a reasonable set of im- 
proved mathematical results, but also to answer certain important questions 
that emerge from such an evidence. 

• The mathematical background for achieving different inequalities related to 
the aforementioned unitary operators has the RS uncertainty principle as 
solid starting-point [12]. The initial algebraic advantage of this investigative 
approach is the inclusion of that contribution coming from the anticommu- 
tation relation between two Hermitian non-commuting operators connected 
via discrete Fourier transform (and/or also related through the Pontryagin 
duality [30]). In fact, this particular contribution allows to include, into the 
algebraic approach, some additional terms — here associated with the mean 
values of different products of the cosine and sine operators — which are 
responsible for correlations between the unitary operators. In principle, the 
theorem derived from this constructive process represents a first important 
point to be strongly emphasized since Eq. (13) introduces a tighter bound 
whose mathematical properties depend on the A^-dimensional state vector 
space in which the initial quantum state is defined. 

• The hierarchical relation (21) derived from the tighter bound characterizes a 
solid bridge between two 'distant' bounds: the first one consists of a zeroth- 
order approximation that confirms the Massar-Spindel result [14] (however, 
it does not depend on the initial quantum state or even on the dimension 
which the state vector space is embedded, these facts being considered as a 
severe limitation for their result); while the second one describes the tightest 
bound for the left-hand side of Eq. (13) and depicts the importance of the 
initial quantum state for a given dimension N. This simple (but important) 
observation justifies our search for finite ground states which quantitatively 
describe those numerical values obtained in Fig. 1 for the tightest bounds. 

• How to construct a finite quantum state which formally explains the tightest 
bounds verified in the numerical calculations? To answer this fundamental 
question, it is necessary to establish an adequate mathematical prescription 
that leads us to obtain a set of normalized eigenf unctions which constitutes 
a complete orthonormal basis in a A^- dimensional state vector space. In this 
sense. Appendix A deals with such a task presenting a reliable mathematical 
relation between Harper functions and tightest bounds by means of ground 
states {|0)Ar} specifically constructed for dimensions A^ G [2, 6]. These finite 
quantum states indeed describe perfectly all the numerical values exhibited 
in Fig. 1 for the tightest bounds, and illustrate the hierarchical relation (21) 
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as well — see Table 1. 
• The bounds determined in this paper for the product Ssu<Ssv exhibit a spe- 
cial link with the theoretical formulation of finite-dimensional discrete phase 
spaces through the discrete Weyl funtion. In fact, this connection establishes 
an interesting link between both the Schwinger and Weyl prescriptions for 
unitary operators, which leads us to guess on the possibility of experimental 
observation with the help of tomographic measurements. 

Now, let us discuss some pertinent points associated with the uncertainty 
principle for Schwinger unitary operators. The first point concerns the Harper 
functions and their remarkable connection with the tightest bounds through 
the ground states {|0)Ar} for a given dimension N fixed. Despite the worked- 
examples in Appendix A belonging to the closed interval 2 < iV < 6, this fact 
does not represent any apparent limitation related to the quantum-algebraic 
framework here exposed. In fact, these results consist of a solid starting point 
for a future search of {{ua\n)}o<n<N-i with a G Z^, whose general expression 
will correspond to a new paradigm for finite quantum states with immediate 
implications in the Fourier analysis on finite groups [19,30] (and/or finite fields 
[31]), as well as in the analysis of signal processing [32,33]. This particular task 
is currently in progress and the results will be presented in elsewhere. 

The second point focus on the systematic study recently developed in Ref. [15] 
for property testing of unitary operators, where D^(\J,Y) := 1 — -^|Tr[U'''V]| 
represents a 'normalized distance measure that refiects the average difference 
between unitary operators'. With respect to this specific measure, Wang shows 
that both the Clifford and orthogonal groups can be efficiently tested through 
algorithms with intrinsic mathematical virtues (namely, query complexities 
independent of the system's size and one-sided error). Since the results here 
obtained describe an uncertainty principle for unitary operators, it seems rea- 
sonable to investigate how these different — but complementary — approaches 
can be juxtaposed in order to produce a unified framework for determined 
tasks in quantum information theory [7]. 

As a final comment, let us briefiy mention that our results also touch on some 
fundamental questions inherent to quantum mechanics (such as spin-squeezing 
and entanglement effects [34]), discrete fractional Fourier transform [35], and 
generalized uncertainty principle into the quantum-gravity context [12]. 
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A Harper's equation, quantum Fourier transform, tightest bound 
and their inherent connections with unitary operators 



Definition (Harper functions). Let {|^)}o<n<Af-i describe a particular set 
of eigenvectors defined in a A^-dimensional state vector space that simultane- 
ously diagonalizes both the Hamiltonian (H) and Fourier {^) operators, that 
is, H|n) = hn\n) and ^\n) = fn\n) for a given dimension A^ fixed. In such a 
case, {hn, fn}o<_n<N-i represents the corresponding set of eigenvalues related 
to the respective Hamiltonian and Fourier operators. Since 

U^\n) = ^U\n) = fnK\n) ^ [H,S^] \n) = 0|n), 

the intrinsic mathematical properties associated with the discrete representa- 
tion {|wa)}o<a<Af_i aHow to obtain the general equation 

(M„|HS^|n) = fnhn{ua\n), (A.l) 

whose solution set {{u^\n)} G M yields a complete orthonormal basis of real 
eigenf unctions genuinely labelled by discrete variables. The Harper functions 
are here attained when one considers H as being a Harper- type operator [32]. 



As a first application, let us, for now, adopt the discrete Fourier operator ' 

N-l 1 Af-1 



/3=0 V N [s^i3'=o 

as well as that Harper Hamiltonian operator previously discussed in Section 
4, namely, H = — sin(6')Cu — cos(^)Cv for ^ G 0, | . Thus, Eq. (A.l) assumes 
the functional form 

N-l 

E €>ia,P;N){u^\n) = /„/i„(M„|n), (A.2) 

/3=0 



where 



0{a,/3;N) = {u^\U\v^) = - 



■ . ,„, f2na\ ,„, /27r/3' 

sin[U) cos [ — — - I + cos(yj cos — — - 



{Ua\Vp) 



represents the mapped expression of the Hamiltonian operator in the discrete 
representations {\ua), I'W/?)} with (-Ualf/j) = -j^u'^^ and u = exp (^); besides. 



^ For A^ odd and discrete labels assuming integer values in the symmetric interval 
[—£, £] with £ = ^^^^ fixed, it is worth stressing that 5^ coincides with that parity 
operator P previously defined in Ref. [12]. 
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{{u^\n)} G M denotes the Harper functions for a given iV G N*. In fact, such 
a result can also be split up into two combined equations as follows: 

/2TTa\ 1 

sin(6') cos ( ^rp j {ua\n) + - cos(6') {{ua-i\n) + {ua+i\n)) = -hn{ua\n) (A. 3) 

and 

^a;"^(M^|n) = /„(M„|n). (A.4) 




The first one describes a three-term recurrence relation and also depicts the 
well-known Harper's equation, whose link with discrete harmonic oscillator 
and discrete fractional Fourier transform was already discussed by Barker and 
coworkers [24]; whilst the second one reflects exactly the eigenvalue problem 
investigated by Mehta [36] when /„ = i" (in this particular case, see Ref. [19] 
for supplementary material), although his ansatz solution 

{Ua\n) = J\fn—7=- Yl ^^P 

V -' ' K= — OO 

does not represent a complete set of orthonormal eigenf unctions [37] — in 
such ansatz solution, A/'n corresponds to a normalization constant and Hn{z) 
denotes the Hermite polynomials. 

Summarizing, Eq. (A. 3) yields, in general, a set of real eigenvalues {hn} whose 
respective eigenfunctions {{ua\n)} constitute a complete orthonormal basis in 
a A^- dimensional state vector space. In this specific case, both the eigenvalues 
and eigenfunctions are dependent on the angle variable 9 G 
us to determine the maximum of cos(^) | (U) | +sin(^) | (V) | ^ 
situation hn ^ {ua\n) with A^ fixed. Thus, the global maximum obtained 
from this mathematical procedure allows not only to fix a given value of 9, 
but also to estimate the smallest eigenvalue of the Hermitian operator H; 
consequently, the corresponding eigenvector will describe the ground state |0) 
characterized by ^ = ^ and /o = +1 for any dimension A^ (it is worth stressing 
that theoretical and numerical calculations confirm these results). Next, let us 
consider the A^ = 2, . . . , 6 cases for 9 = j fixed, in order to provide a complete 
list of results exhibited in Table 1 associated with the normalized ground state. 

• Case N = 2 (prime dimension) and ho = —1. This first example obeys the 
criterion "easy to calculate" , once the corresponding ground state 



0, 1 , which leads 
br each particular 



/2 + v/2 , , , V2-v/2 , 
10/2 = - — t; po) + — n pi) 



7 



Such a quantity implicitly defines the boundary — or, more precisely, the convex 



hull — of the accessible region related to the {|(U)|, |(V)|}-space [14]. 
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allows to obtain '^(v) = | and "^uC^v) = 1, which implies that S^J = 1. 
Case A^ = 3 (prime dimension) and ho = — v6+^^ . This particular example 
yields certain interesting peculiarities since its respective ground state 



|0)3 = y— ^l^o) + V^2~'''^^ + V ^2~'"^^ 

leads to achieve '^(v) = | + ^~g and 1^u((5v) = 15 — 8v^, which corroborate 
those results obtained by Opatrny [23] for the number and phase operators. 
Furthermore, it is easy to demonstrate that S^J = 7 — 3\/3, which justifies 
the approximated numerical value 9^i ~ 3.254 exhibited in Table 1. 
Case N = 4 (even dimension) and Hq = —1. In such a situation, the ground 
state 

10)4 = — ^ — \uo) + — Fi) + — ^ — F2) + — Fs) 

attains those same values of "^(v) and %ij(sv) verified in A^ = 2; however, it 
is worth stressing that S^J = 2, since e = | in this particular case. 

non- 



Case N = 5 (prime dimension) and ho = — — ^j^^ — ^'^. For this 

trivial example, the normalized ground state assumes the form 



10)5 = , — (|mo) + dil^i) + 021^2) + 02!%) + ai|M4)) 

"+2(a? + ai) 



with Oi = ^^— ^^ — and 02 = ^^-^ — . So, after some lengthy 

calculations, the exact expression "PsviSv) = — ~ ~^^ for the variance 

19+v^+V110+38v^ 

allows to show that 



5 + V25 + 10^5 j f 45 - V5 - V 110 + 38^5 
5 f 19 + V5 + \/llO + 38V5 



S^'J = ^ , '-^ , '- ^ 1.9444, 



corroborating, in this way, the numerical result 9^i ~ 3.781 (see Table 1). 



Case N = Q (even dimension) and ho = —— — ^ . Similarly to the previ 

ous case, in this situation the ground state admits 



|0)6 = , i\uo) + bi|Mi) + b2|M2) + bslws) + b2|M4) + biJMs)) 

^i+2(bf + bi) + bi 



as a solution of Eq. (A. 2), where 



27 



-2 + V5 + V2T -A + VQ + VU 

bi = = 

2 4 



5 + 721 - 3V5 + ^21 10 - SVe - 3y/U + 2y/21 

bo = = 

^ 2 4 



4 + Ve - 721 + 2V5 + 721 = -4 + 2Vq + v/l4 - ^21 



represent the respective coefficients. Note that ysv{5v) = 19 — 4v21, which 

imphes in S^^^ = 19(1 + ^3) - 4^42(2 + x/3) ^ 1.8297 (this resuh justifies 
that numerical value 9^i ^ 3.348 appeared in Table 1 for N = 6). 

As an initial purpose, these first theoretical results related to the ground state 
are sufficient to clarify the numerical results depicted in Fig. 1 for the tightest 
bound. In fact, the results exhibited in this appendix indeed represent a first 
investigative step towards a general mathematical recipe that presents as a 
primary product the eigenfunctions {('Ua|^)}o<ra<Af-i, which differ from that 
Mehta's ansatz solution. 
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